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Abstract 

We compute the spectrum of /3-decay, assuming that (Majorana or Dirac) neutrinos prop- 
agate in constant potentials. We study the modifications of the spectrum due to the effect 
of these potentials. Data on tritium decay and on '^H-^He mass difference allow us to infer 
bounds in the electronvolts range on the potentials. 



We study the modifications of tlie electron spectrum in /?-decay due to liypotlietical 
potentials influencing the electron neutrino propagation. The weak interactions of the 
neutrinos with ordinary matter do modify the neutrino propagation j^, but are not 
expected to be able to affect the measured /3-spectrumQ. In the present work we take 
a purely phenomenological point of view, namely we do not attempt a study of the 
nature of the potentials, but just wonder which kind of observable effects they could 
produce. Apart from the interest in contemplating a theoretical possibility, some 
additional motivation for this work comes from anomalies suggested by experimental 
studies of tritium decay [0, Q . 

We show that potentials in the electronvolt range can lead to observable features in 
the /3-spectrum, discuss the features of the spectrum near and far from the endpoint, 
and use the experimental data to put an upper bound on these potentials. 

The structure of the paper is the following: in the first section we review the 
formalism for treating Majorana or Dirac neutrinos (essentially equivalent to the 
formalism put forward in for Dirac neutrinos); in the second section we compute 
the /3-spectrum in the presence of a background potential for the neutrino; in the last 
section we discuss the features of the spectrum, and compare our results with the 
experimental findings and with other theoretical models. 



1 Neutrino fields in presence of background po- 
tentials 

Consider a neutrino described by the bi-spinor A^^ = 1,2), that participates in 
weak interactions and propagates in a constant potential Vl : 

C% = iNaad'^N -{mNN + h.c.) - Vl NN; (1) 

where a = 0, 1, 2, 3 are Lorentz indices and m is a Majorana mass. The spinorial 
contraction^ imply that all terms but the last one are Lorentz invariants; therefore 
this latter term implies the existence of a preferential Lorentz frame. Alternatively, 
assuming that a second bispinor A^^ exists, we can consider the propagation described 
by the lagrangian: 

£0, = iNoad^N + iWaad^N" - (m A^^A^ + h.c.) -VlNN- VrWN", (2) 

where m now denotes a Dirac mass term. A^^ will be assumed not to participate in 
the weak interactions. Eq. (|^) is invariant under those U(l) transformations that act 

*A density of scatterers p of the order of the Avogadro number gives origin to a potential Vweak 
Gpp ^ 5 • 10^^'* eV, 14 order of magnitudes smaUer than the present experimental sensitivity; 
compare however with |^ . 

^^Our conventions for bispinors: Na — N^; N"^ — c^^Nb (e is asymmetric, and e^^ — 1); 
NN = N'^Na; Oa = (1, o) [a are the Pauh matrices); Noa^ = N^a^^Ne; NN = Nj^Na- 
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with opposite charges on and A^*^ (Majorana mass terms hke NN and/or potential 
terms hke A^'^A^ exphcitly break this invariance). Introducing the four-spinors: 



^f=(fl) [Majorana] and = ) [Dirac], (3) 

we can rewrite the lagrangians (|l|) and (Q), after performing a partial integration in 
the corresponding actions, in the following form: 

^M = l^ i^lad'' -m + 7%Vl] u^' and 

2_ (4) 

= [ilad'' - m - 7°(r5 - 7^^)] z/^. 

We used the gamma matrices in the chiral representation (but of course the la- 
grangians have the same form in any representation), and set the definitions: 

Vl + Vr 



s — ^ 

Vl-Vr 

The solutions of the equations of motion are easily obtained. In the Majorana 
case we search for the positive energy solutions of the form: 

u^\x) = exp[—iE{X)t + ip-x] x u{X,p); (6) 

where A = ±1 denotes the helicity: T, ■ p / p [u{X, p)] = X u{X,p), using the notation p 
for the squared three-momentum. Since the spinor u obeys: 

[E(A)7° - m - 7V(Ap - Vl)] u{X, p) = (7) 

the potential can be taken into account replacing 

P^P-XVl (8) 

in the usual formulae for energy and for spinors. In particular, the dispersion relation, 
depicted in fig. 1, becomes: 



E{X) = ^w? + {Xp - VlY. (9) 

The steps (^-|^) can be repeated for the Dirac field, since after the simple replacement: 

Vl Vd (10) 

the wave v'{x) = exp(iVst) u^{x) is shown to obey the same equation of u^\x). 
The quantization of the Majorana field yields us: 

E / k^'^ "(^'^) e-^^W*+*^^-"^ (11) 



p'Ux) 



2E{X){2n) 



+ia\X,p) Cu{X,p) e 



+iE{X)t—ip-x 
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where spinors and operators are normahzed to 2E{X), and ^ is a phase factor entering 
the Majorana condition: Cu^^i^x) = ^v^^x). In order to obtain the quantized Dirac 
neutrino one has {%) to replace the creation operator in (|ll|) with a different 
degree of freedom h\ (ii) perform the substitution (|10|) into the expressions for en- 
ergy and for spinors, and finally {in) multiply the resulting field u' by the phase 
exp{—iVst) . Due to this last step, Vs will enter in the expression for the energy of a 
Dirac neutrino (resp. antineutrino) as: E{X) + Vs (resp. E{X) — Vs). 



2 /3-spectrum 

The /3-decay is due to the charged current lagrangian: 

Ccc = ^^'^(l - ^^)d e7a(l - (12) 

where = z/^ or u^. The potentials modify the /3-spectrum. In this section we 
describe in detail the modifications which take place for a Majorana neutrino, and 
indicate how to adapt the formulae to the Dirac case. To be concrete, we will focus 
on the tritium decay: 

^H^^He + e + I/, Q = 18590 eV (ref. [6]) (13) 

for which the energy release Q (the nuclear mass difference being AM = rrif. + Q) 
is so small that all particles except possibly the antineutrino can be considered non- 
relativistic. 

The matrix element squared is: 

\M{\)\^ = const. X {E,{\) + \p, - Vl), (14) 



where the constant depends on the fundamental parameters in ([T^), on the nuclear 
matrix element and on the masses of the electron and of the nuclei (one simple deriva- 
tion of this result: consider the usual result for |A^p for neutron decay, substitute 
(py)" with (pi/)" + Xrriy (n^)" to account for antineutrino polarization, perform the re- 
placement d^) and then take the non-relativistic limit). For the Dirac case, use (p!0[). 
As expected, the matrix element corresponding to the positive helicity dominates in 
the relativistic limit. 

Since the antineutrino energy depends on the polarization, the maximal electron 
energy will also depend on it. The minimal antineutrino energy is E'™"' = my if 



A = sign(VL), or i?™'" = yml + Vl if A = — sign(VL). Therefore, from the energy 
conservation condition we obtain: 

r AM -m^ if A = sign(yL) 

Er''W = { + vl if A = -sign(Vz.). ^^^^ 
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The first case corresponds to a non-zero minimum momentum of tlie antineutrino: 
Pu = 1^1- For tlie Dirac neutrino, in order to obtain tlie maximum electron energy 
one lias (i) to perform the replacement (|1^) and (ii) to add Vs in the previous formula. 

Let us consider the phase space of the reaction. Using the delta function from 
momentum conservation to integrate away the helium-3 three-momentum, integrating 
over the possible direction of the electron, and averaging over those of the antineutrino 
we get: 

where F^{X) is a phase space factor related to the antineutrino degrees of freedom. 
When integrating over the unobservable antineutrino momentum, one has to take 
into account the dependence of the energy on the helicity: 

(1) If Ee < AM — J ml + Vl (corresponding to > J ml + V^, see fig. 1), there 



are two possible momenta, one for each helicity state: = E^ — ml + XVl- For 
them we get: 



-ml + XVl^^ 



E?. — m?, 



(17) 

=AM-E. 



(2) In the case E^ > AM — J ml + Vl there is only one helicity state, A = sign(V/^) 



but there is a second value of the momentum at fixed energy: py = El — ml 
+\Vl. The phase space factor for this determination of the momentum is: 



Fu{\) 



El - ml - XVl 



\2 



E?, — m?. 



(18) 



whereas, for the other momentum, formula (^Tjj applies (with the same value of the 
helicity). Notice that this phase space factor is formally identical to the one that we 



would obtain using the A = — sign(Vi) and Pi, = ^ El — ml + XVl, even if, strictly 
speaking, this "square momentum" p^ is negative. 

The differential spectrum of the electron can be estimated by summing up the 
antineutrino helicities, dT = J2x d^si^X) |A^(A)p / (2 msHe), where as discussed above, 
both the matrix elements (0) and the phase space (|16D depend on the helicity: 



^p^E^F{Z,E^) X + {Vy + aVLr 



dE, — ^ ' - 2V, 



(19) 

^AM-Ee 



in this formula p^ = E^ — ml and Vy = ^ El — ml; F{Z, Eg) embodies the electro- 
magnetic corrections due to the interaction of the final states. The maximal electron 
energy is 

^ondp. ^ _ [Majorana], (20) 
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and is not shifted by the potential Vl (we recall that AM denotes the '^H-^He mass 
difference) . 

In the Dirac case we get: 



^ oc p^E^FiZ, E^) X E J'^^^^" + a Vn f 



(21) 

E^=AM-Ee-Vs 



The maximal electron energy is 

^ondp. ^ _m^^Vs [Dirac], (22) 

and therefore differs from the Majorana case ii Vs ^ 0. 
The Kurie plot is defined as: 



if masses and potentials are absent, K{E) turns out to be a line approaching zero at 
the endpoint. In the case of small neutrino mass the potential enters the formula for 
the Kurie plot as: 

K{Ee) oc I AM -Ee + Vil, small (24) 

for both Majorana and Dirac cases, but the endpoint is in a different position unless 
Vs = (compare eqs. (|0D and (|2|)). 

The integral spectrum, often used in the experimental analysis is: 



3 Discussion of the results 

We first discuss some typical features of the /5-spectrum, and then compare them 
with the experimental results and with other theoretical models. 



3.1 Features of the /3-spectrum 

We start with some general remarks on the differential spectrum, eqs. (0) and (pT]): 

1. While the effects of the masses on the spectrum are of order {m^/Eiy)'^, the 
effects of the potential are of order V/E^. 

2. The Majorana and the Dirac neutrino spectra for equal masses are identical, 
not only if the potentials are absent, but under the weaker assumption Vs = 
(which, by definition (^, is equivalent to Vd = Vl). 
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3. There is a curious behaviour of the differential spectrum at the endpoint: If 
the mass and the potential are both non-zero, it goes to infinity, since — > 
in eqs. (|19|) and (0). This fact however does not imply contradictions with 
the physical interpretation, since in any given finite interval the probability of 
transition is finite. This behaviour is due to the increasing number of states 
corresponding to a fixed energy interval when the neutrino energy decreases 



toward its minimum, compare with eqs. (|T6|-0), or see fig. 1 (we recall that 
the minimum corresponds to the-infinite-set of states = \Vl \ n^,, where the 
orientation is arbitrary)]^ 

Continuing the discussion of the differential spectrum, we come to the differences 
between the Majorana and the Dirac cases. For a Majorana neutrino it is sufficient 
to distinguish if Vl is larger or smaller than zero. For small mass, in the first case the 
spectrum just translates upward, whereas in the second case it translates downward 
for energies < AM + Vl, where it becomes zero, and then it increases up to the 
endpoint (pO|): a dip is present. These features are clearly illustrated by eq. (^). The 
Dirac case is identical, apart for the fact that the position of the endpoint depends 
on Vr. Notice in particular that for Vj, = Vr (= Vs, Vd = due to eq. (^)) the Kurie 
plot is simply a straight line which approaches zero at = AM + Vl- 

Let us conclude with some comments on the integral spectrum. In the Majorana 
case, the integral spectrum for Vl > moves to the left, and approaches in a steeper 
fashion the endpoint when the mass is increased. Instead if Vl < 0, for masses of 
the order of \Vl\ or smaller the integrated spectrum remains approximatively con- 
stant in the region that extends 2\Vl\ below the endpoint (the size of the dip in the 
differential spectrum). In the Dirac case for positive and increasing values of Vr the 
endpoint moves to the right, extending the region in which the integral spectrum 
stays approximatively constant, in accord with the features already discussed. 



3.2 Comparison with the experimental results 

Experimentally one can study the shape of the /3-spectrum in the region far from the 
endpoint, or near the endpoint. Let us discuss these two regions in turn. 



3.2.1 Region far from the endpoint 

Due to the first feature of the spectrum emphasized in the previous subsection, a mass 
of some eV would give a deviation of the spectrum only near the endpoint, whereas a 
potential Vl of few eV would give a detectable shift of the spectrum also in the region 
far from the endpoint (the eV range being the present experimental accuracy). In this 

^li the potential Vl is zero, the minimum shifts at p^, = 0, and therefore the phase space factor 
dpi, shrinks to zero when the minimum is approached, whereas if = the dispersion relation 
curve in this limit docs not flatten, so that the phase space attains a finite value. 
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hypothesis, the experiments which measure the ^H-^He mass difference fitting the j3- 
spectrum far from the endpoint are sensitive to AM|^ = AM + Vl (eq. applies 
since the neutrino energy is much larger that its mass in the region under discussion) 
which would systematically deviate from the AM|j^q^_^, obtained by techniques of 
measurement of the mass difference which do not rely on the /5-spectrum. Whereas 
older data were to some extent indicative of a discrepancy between the /3- and the non- 
/3-determinations of AM, more recent and precise determinations are in remarkable 
agreement 0. Including the six more recent measurements reported in ref. [0] (5 (3- 
and 1 non-/5-techniques) we estimate: 

Vl = AM|^ - AM|^„„_^ = 0.3 ± 2.7 eV (26) 

Eq. (HI) rules out Vl > 5.6 eV or Vl < -5.0 eV at 95% CL. 



3.2.2 Endpoint region 

The region near the endpoint is very interesting for the search of electron antineutrino 
mass. Fits to recent data ||, || suggest a "negative masses squared", ml ~ —20 
eV^. Recalling the well-known feature that a (positive) neutrino mass term depletes 
the endpoint region of the differential spectrum, one realizes that a "negative mass 
squared" simply amounts to the fact that the measured spectrum is even larger than 
what one would expect postulating a massless neutrino. Putting aside the possibility 
of statistical or instrumental effects, such findings could be taken as suggestions of a 
peak-like structure near the endpoint region. 

In the following we focus on the possibility to use the model under discussion to 
explain this feature. But before doing so, let us spend a few words of warning: 

• Even if certain choices of the potentials are able to explain the presence of a 
peak in the endpoint spectrum, one should keep in mind that different types of 
patterns can be obtained with other choices of potentials (and, in particular, of 
their signs). In addition, other models exist which may entail modifications of 
the endpoint spectrum (see also the discussion below). 

• The fit value: ^ —20 eV^ quoted above is dominated by the experimental 
findings of the Troitsk collaboration [Q, which reports the smallest errors, 4.8 
and 5.8 eV^. The confidence level of the average fit values of and @| is of 0.6% 
or 20%, if we sum quadratically or linearly the statistical and the systematic 
errors (when they are quoted separately). These considerations suggest that it 
will be quite important to follow future experimental developments. 



As it is clear from eq. (24), a Majorana neutrino with a small mass and a negative 
potential Vl of few eV could explain the appearance of a peak in the differential spec- 
trum exactly at the endpoint, and a dip \ Vl\ electronvolts below. But for a Majorana 
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neutrino the position of the dip is hnked to a downward shift of the spectrum in the 
region far from the endpoint, and therefore the bound ( P^D apphes. 

This hmitation does not hold necessarily for a Dirac neutrino. In fact, if V/, ~ 
and Vr is positive, there is no deviation in the differential spectrum up to the point 
in which it vanishes, Eg = AM, but for larger energies the spectrum would start to 
rise constantly, up to the endpoint Eg = AM + Vr/2 (we assumed again a small 
neutrino mass term). We conclude that a potential Vn of some eV, felt by a Dirac 
neutrino, can manifest itself in a peak of the differential spectrum near the endpoint 
region. A Vr in the twenty eV range or so can already be excluded by the data. Let 
us stress that in such a scheme there is no need to attach any physical meaning to 
the "negative mass squared". 

3.3 Comparison with other theoretical models 

An observable deviation of the spectrum in the region far from the endpoint region 
(eq. ( |26D and related discussion), characterizes the scenario in which the neutrino 
propagates in a (large) background potential Vl- 

Instead, peak-like structures in the endpoint region can also arise in the presence 
of a degenerate neutrino sea, filled up to the Fermi level /iF [H- In fact, the reaction 
of absorption of a neutrino z/gea + ^H e + '^He permits the emission of electrons with 
energies between AM and AM + /i_F, if the neutrino mass is supposed to be small 
in comparison with ^p. Therefore, one cannot discriminate between these two theo- 
retical possibilities from the experimental information on the /3-spectrum. Instead, 
if neutrino masses are not small it is possible, at least in principle, to distinguish be- 
tween the two cases {i.e. potential versus chemical potential). In fact, if we suppose 
the existence of a degenerate sea of massive neutrinos, no electron could be emitted in 
the energy range between AM — rrty and AM + rriy (the emission due to the neutrino 
absorption taking place for larger energies, up to E^. = AM + rriy + fip) 0. Such 
a discontinuity of the differential spectrum is absent in the models discussed in the 
present work. 
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